% soit le systeme linéaire 
% 1x1+2x2+3x3 -4x4=-2
% 5x1+6x2+2x3+x4=27
%3x1+4x2+5x3-3x4=14
% 2x1+7x2+2x3+x4=26

a=[1 2 3 4;
   5 6 2 1;
   3 4 5 3;
   2 7 2 1];
X=[-2 27 14 26];
b=[a*X];
n=size(a,1);
for k=1:n-1
{
for i=k+1 : n+1
{
w=b(i,k)/b(k,k);
for j=k : n+1
b(i,j)=b(i,j)-w*b(k,j);
fprintf("%d",b(i,j));
}}
for i=n-1 : 1
{
s=0;
for j=i+1 : n
{s=s+b(i,j)*X(j);}
X(i)=(b(i,n+1)-s)/b(i,i);
fprintf("%f" ,x(i));
} 